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ABSTRACT 

ESO 566-24 is an extraordinary barred galaxy, which has four regularly spaced spi- 
ral arms in blue light images. This type of four-armed spiral structure, which is rare 
among the spiral population, is clearly seen also in near-infrared images, and thus 
is present in the old stellar population. We have constructed dynamical models of 
ESO 566-24 in order to understand the cause of its four-armed structure. The disk 
gravitational potential is determined from near-infrared photometry, and the gas dy- 
namics is modelled as inelastically colliding particles. The resulting morphology and 
kinematics with different assumed pattern speeds, disk vertical thicknesses and dark 
halo contributions is compared with observations. Our models can reproduce the main 
morphological features of this galaxy: the four-armed spiral, and the inner and nuclear 
rings. The pattern speed of the bar in this galaxy is such that the corotation resonance 
is well outside the bar radius, rcti/Har = 1-6 ± 0.3. The four-armed spiral resides in 
the region between inner and outer 4/1-resonances. Also, the main kinematical fea- 
tures, including bar-induced deviations from circular rotation, are explained by our 
models. The best fit is obtained when the dark halo contribution is just enough to 
make the modelled rotation curve match the observed one. This "minimum halo" is 
rather moderate: luminous matter dominates the rotation curve within the disk region. 

Key words: galaxies: evolution - galaxies: fundamental parameters - galaxies: kine- 
matics and dynamics - galaxies: spiral - galaxies: structure - galaxies: individual (ESO 
566-24) 



1 INTRODUCTION 

Barred galaxies constitute the majority of all spiral galax- 
ies. According to recent near-infrared (NIR) studies, about 
70% of disk galax ies have a bar ( Es kridee et alJ l2000t 
iKnapen et all [2000 ) . Barred galaxies often have rings or 
pseud orings in their disks w ith very distinctive character- 
istics feuta fc Combeslll996l) . Depending on its relative size 
with respect to the bar, a ring can be classified as an outer 
ring (semi- major axis about twice the bar radius), an inner 
ring (semi-major axis w bar radius) or a nuclear ring (semi- 
major axi s about 10 % of the ba r radius, though the scatter 
is large) (iButa fc Crockeilll993l) . Although all three types 
of barred galaxy rings can show enhan ced star formation 
(see, for example, IButa fc Purcell||l998l), nuclear rings are 
most prone to starbursts iMaoz^tallT996l : lButa fc Combes! 
2000). In addition to rings, barred spiral galaxies often show 
offset dust lanes on the leading sides of the bar (assuming 
trailing spiral arms). 

* E-mail: pertti.rautiainen@oulu.fi 



The dynamical behaviour of barred galaxies has been 
modelled in different ways. Some studies are general, in the 
sense that no specific barred galaxy is modelled, but they 
concentrate instead on issues such as the formation and evo- 
lution of the bar. In most of these studies, the large-scale 
bars are suggested to be fast rotators, i.e. their corotation 
resonance (CR) radius should be near the end of the bar 
<SeUwoodll98lfl . This is in accordance with gas dynamical 
simulations of the shap es of the le ading offset dust lanes 
inside the bars llAthanassoulall l992bl). E xceptions to fast ro - 
tati ng bars are the studies bv lCombes fc Elmeereenl dl993fi 
and lMiwa fc Noeuchil il998fl ; in the former study a short and 
slow rotating bar forms in a simulation where the disk was 
strongly truncated radially, while in the latter, a slow bar 
forms as a result of an interaction. The contribution of a dark 
halo inside the radius of the optical dis k can have a strongj 
effect on disk dynamics. According to lOstriker fc Peebles! 
l|l973|) a massive halo can inhibit the bar formation, and 
iDebattista fc Sellwoodl |l998, 2000|) have shown simulations 
where the interaction between a bar and a halo decelerates 
strongly the bar rotation. However, conflicting results have 
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also been published : lAthanassoula fc Misiriotisl i2002T) and 
lAthanassoulal (2002) show exa mples where the halo a c tually 
promotes bar formation, and IValenzuela fc Klvpinl (2002) 
found a rather low angular momentum transfer rate from 
the bar to the dark halo. 

In principle, the pattern speed of a bar can be de- 
termined from observations for a disk wit h a well-defined 
pattern speed jTremaine fc Weinberdll984h . The Tremaine- 
Weinberg method further assumes that the used tracer com- 
ponent follows the continuity equation, i.e. its intensity is 
proportional to total density. When these conditions are 
fulfilled, the pattern speed can be determined from the 
observed intensity and the line-of-sight velocity along a 
strip parallel to galaxy major axis. In most of the cases, 
where this method has been applied, the cor otation res- 
onance seems to be near the end of the bar (iKentJ fl987t 
iMerrifield fc Kuiikerll995tlAguerri et alJliool) . Note, how- 
ever, that due to its limitations, this method has been suc- 
cessfully applied mostly to SBO-galaxies. 

Another kinematic al method to d etermine the pattern 
speed was presented bv lCanzianl il993l) and is based on ob- 
served residuals in the velocity field after subtracting the ax- 
isymmetric velocity component. In the case of a two-armed 
spiral, the velocity residuals exhibit an m — 1 spiral shape 
inside corotation, whereas outside corotation, an m = 3 spi- 
ral should appear. This method was used to de t ermin e the 
pattern speed in NGC 4321 by ISempere et alJ il995T) and 
ICanzian fc Allenl il997t) . who found a rather low value, set- 
ting the corotation radius about 1.6 - 1.8 times the bar ra- 
dius. However, the method is very sensitive to noise, and the 
construction of the axisymmetric velocity component is not 
trivial. Purcell (1998), who studied a larger galaxy sample, 
found that only in 7 of 27 Sbc galaxies could the expected 
pattern be seen in velocity field residuals. 

There are several suggested resonance indicators related 
to morphological details of spiral structure dElmegreen et all 
I1989T) . Even if these are valid, their application to barred 
galaxies is not necessarily straightforward: there can be 
more than one mode in the system. For example the outer 
spiral arms can have a l ower pattern speed than the bar 
JSellwood fc Sparkdll98l) . In such cases the pattern speed 
derived e.g. from the dust lane morphology of the outer arms 
would not correspond to that of the bar. 



The formation of rings in barred galaxies has been 
explained by the response of the dissipati ve gas compo- 
nent to the gravitational torque of the b ar ( Schwarz 
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Bvrd et alJ 11994 iButa & Combes! 
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2000). The shapes of the rings re- 



semble the shapes of certain closed period ic orbits in barred 
potentials JContopoulos fc Grosb0lll98gb . When the non- 
axisymmetric perturbation is weak, these rings form very 
close to resonance radii. The usual resonance identifications 
are such that the outer rings are located near the outer 
Lindblad resonance (OLR), the inner rings near the inner 
4/1-resonance and the nuclear rings near the inner Lind- 
blad resonance (ILR). Thus if the rotation curve of the 
galaxy is known, the sizes of the rings can be used to es- 
timate the bar pattern speed. However, these pattern speed 
estimates can be misleading when the bar perturbation is 
strong, and the resonance distances based on the linear ap- 
proximation fail. This is especially the case with nuclear 
rings, which are located in a region where the relative non- 



axisymmetric forcing is strong. Furthermore, both nuclear 
and outer rings can often be affected by another mode in 
addition to the main bar; nuclear rings by a secondary bar, 
and outer rings by an outer s piral mode IIRautiainen fc Said 
l200dlRautiainen et alJl2002h . 

The effect of the strength of the bar is taken into 
account when detailed models for individual galaxies 
are being constructed. The adopted approach has var- 
ied fr om fitting an analytical ma ss model to observa- 
tions jDuval fc Athanassould I1983T) to d etermining the 
mass model from near-infrared o bservations llLindblad et alJ 
Il99d iLindblad fc Krister]|l996h . In some models, the non- 
axisymmetric forcing of the bar was not enough to repro- 
duce the spiral structure with the observed ext ent, and 
so an unseen ova l-shaped component was added jEngl and 
119891 iBalll I1992T) . A spiral potential corresponding to 
the NIR-m orphology solved this problem more elegantly 
llLindblad et al]ll996l) . The number of successfully modelled 
barred galaxies is still quite small, and it is not clear how 
well the ranges of parameters such as the bar pattern speed 
and the relative strength of the non-axisymmetric perturba- 
tion are known. 

We started our modelling project of ringed galaxies 
with the normal, but r ather weakly barr ed galaxy IC 4214. 
From the observations (IButa et alJll999l ), we constructed a 
mass model with a few free parameters JSalo et al.lll999Fl . 
In our best-fitting models both the morphology and kine- 
matics were well reproduced. Following the successful mod- 
elling of IC 4214, we extend our studies here by constructing 
a dynamical model for an extraordinary barred spiral: the 
four-armed galaxy ESO 566-24. 



2 A REMARKABLE FOUR-ARMED BARRED 
SPIRAL 

ESO 566-24 is a barred spiral galaxy of de Vaucouleurs re- 
vised Hubble type SB(r)b. It was first recognized as being 
an exceptional example of a four-armed ringed barred spiral 
galaxy during early stages in the produ ction of the Cata- 
log of Southern Ringed Galaxies (CSRG, lButalll995h . High 
quality optical CCD images in the U, B, and I passbands, 
obtained in follow-up work bv IButa fc Crocked (Il99lf) . re- 
veal an extremely regular four-armed pattern of normal spi- 
ral arms, as well as star formation in the inner ring and 
in a very small previously unknown nuclear ring. These im- 
ages also revealed a faint, partial outer rin g-like pattern tha t 
envelopes the four-armed spiral pattern jButa et al.lll998T) . 
Although rings are not rare in barred galaxies, regular four- 
armed outer spirals are very rare and, if seen at all, appear 
in a partial manner, as the second ary spiral a r cs and two 
main outer arms seen in NGC 1433 iButal ll986: But a et alJ 
2001). 

IButa et alJ (Il99gf) showed that the four-armed struc- 
ture of ESO 566-24 is clearly seen in a NIR //-band image. 
Thus, the pattern is also present in the old stellar popula- 
tion. A r ecent study of about 2 00 galaxies in different wave- 
bands by Eskr idee et alJ a2002i) demonstrated that the NIR 
morphology of galaxies is usually not as dramatically dif- 
ferent from the optical morphology as is sometimes stated 
(Bloc k~fc Puerarl I1999T) . with multiple-arm and flocculent 
morphology also being seen in the //-band. However, in most 
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of the cases, two arms dominate, and the extra arms are 
weaker and usually form as a bifurcation of the main arms. 
This is not the case with ESO 566-24: the four arms are of 
about same strength and length. 

In this paper, we use numerical models to find the pa- 
rameter range (pattern speed, mass-to-light ratio) which can 
reproduce the regular four-armed morphology of ESO 566- 
24. We study how resonances are located with respect to 
the observed structures such as rings, and we also identify 
orbits supporting these structures. Finally, we compare the 
spectacular morphology of ESO 566-24 with more partial 
patterns seen in other barred galaxies. 



3 METHODS 



The method used by ISalo et, al] (|l999F) in modelling the 
early- type, weakly-barred spiral galaxy IC 4214 is used here 
for modelling ESO 566—24. The gas morphology is mod- 
elled by inelastically colliding test particles, which move in 
a rigidly-rotating gravitational potential. We estimate the 
mass distribution from NIR ff-band photometry, assuming 
a constant mass-to-light ratio throughout the disk. Before 
deprojecting the galaxy image, a bulge-disk decomposition 
is done, and the resulting bulge component is subtracted 
from the image. Then the disk surface brightness (and sur- 
face density) distribution is approximated using a Fourier 
decomposition 



E(r,0) = E (r) |l + ^ 4 m (r) cos [rri (0 - m )] j 

OO 

= E o (r) + ^E m (r,0), (1) 

m — 1 

where r and 9 are the polar coordinates, E(r, 9) the disk 
surface density, Eo the axisymmetric surface density, and 
A m (r) and # m (r) are the Fourier amplitude and phase an- 
gles. For our analysis, we used even components m = - 
8. This introduces smoothing of small scale structure such 
as star forming regions but has little or no impact on the 
basic backg round stellar potent ial. A similar procedure was 
adopted bv iKranz et alJ a2001f) . From the density compo- 
nents E m , we then calculate the potential components "!> m 



$ m (r,(9,z = 0) 



-G 



r'dr' / E m (r',e')g(Ar)d9', (2) 



where (Ar) 2 = r' 2 + r 2 — 2rr'cos(9' — 9) and the convolution 
function 



g(Ar 



gz (z)((Ar) 2 + z 2 )- 1/2 



dz 



(3) 



is used for calculating the effect of the vertical density profile 
q(z) (normali zed to unity when integrated over the vertical 
direction). In ISalo et all jl999F) the gravity softening was 
used as a measure of finite disk thickness. In this study we 
use an exponential vertical distribution, 



Q(z) = ^^p{-\z\/h z ) 



(4) 



which has a similar effect. A more de tailed description of 
the potential calculation can be found in lLaurikainen fc Said 



(2002). We have tried scale heights in the range 1/12 - 1/2 
times the radial exponential scale length h r of the disk, cor- 
responding to the change from a very thin disk to an unre- 
alistically thick one; for the m orphological type of ESO 566- 
24, h z /h r » 1/6 is expected Jde Griis et al.lll997l : Tde Griisl 
1998). Some tests were made with variable scale height, mak- 
ing the disk thicker in the bar region. 

The non-axisymmetric perturbation of the potential is 
turned on gradually during four bar rotation periods. This 
procedure was adopted to smoothly "guide" particles to re- 
alistic orbits, thus trying to minimize the effect of transient 
phenomena related to abrupt changes in the potential, e.g. 
such as a large number of gas particles ending up on crossing 
orbits. 

The gas component is modelled as a two-dimensional 
disk consisting of 20,000 inelastically colliding test parti- 
cles. In each collision, the relative velocity component par- 
allel with the line joining the particle centres is reversed and 
multiplied by the coefficient of restitution. In each time step, 
only one collision per particle is allowed. The particle mo- 
tion is integrated using leap-frog integration. The time-step 
is chosen so that at the radius of the smallest morphological 
detail, the nuclear ring, one circular rotation period is about 
50 steps. In addition to gas particles, we have also made sim- 
ulations with 200,000 non-colliding "stellar" particles with 
different initial velo city dispersions. F or more details on the 
simulation code, see lSalo et al.l (|1999). 

We use a modifi ed version of the bulge model presented 
in IButa et alJ lll998f) . We take the average radial profile of 
the bulge, assume the bulge to be axisymmetric and calcu- 
lat e its g r avitat ional potential using the algorithm described 
in IKend Il986'l. As a possible halo component we use an 
isothermal sphere with a constant core radius, which has a 
rotation curve 



Vh(r) = v a 



(•>) 



where Vca is the asymptotic velocity at infinity and r c is the 
core radius. 



4 OBSERVED CHARACTERISTICS OF ESO 
566-24 

B ecause the observa tions and data reduction are described 
m IButa et al.Hl998ft . we here just give the main morpholog- 
ical and kinematical characteristics of ESO 566-24, concen- 
trating on details which are useful for comparing simulations 
with observations. 

ESO 566-24 has a clear bar component in both visual 
and NIR-images. The bar is considerably stronger than in IC 
4214. It is difficult to judge what would be the best method 
to determine the bar length. Applying several methods to 
_ff-band images: e.g., axis ratios of isophotes and local am- 
plitude minima of m = 2 and m = 4 Fourier components 
of disk surface brightness, we estimate the bar radius to 
be ri, ar = 18" ± 2", or about 1.6 times the exponential 
scale length of the disk, h r ps ll'.'3. T his is in rea s onabl e 
agreement wi th the value adopted in | B uta et alJ l)l99S^ . 
rbar = 16'.' 7. lAthanassoula fc Misiriotisl l|2002T) use seven 
different methods to determine the bar length in N-body 
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Figure 1. B- and H-band images of ESO 566-24. The images on the left side show the galaxy as it is oriented on the sky (foreground 
stars have been removed) and the images on the right side have been first deprojected and then rotated so that the bar is horizontal. 
The width of the frames is 120". 



simulations. However, we found that the complicated mor- 
phology of ESO 566-24 made many of these unsuitable (e.g. 
their m — 2 phase angle criterion would give the bar radius 
in the middle of the spiral pattern). Colour index maps, e.g. 
B — I, sho w that there are dust lanes on the leading sides 
of the bar jButa fc CrockeJll99ll) . 

The bar is surrounded by an inner ring, which appears 
circular in the plane of the sky but deprojects to an ellipse 
parallel with the bar major axis for the adopted orienta- 
tion parameters: position angle PA — 73° and inclination 



i = 43° iButa et al.lll99gl) . Parallel align ment is a n ormal 
characteristic of barred galaxy inner rings feutalll995l) . The 
deprojected semimajor axis radius of the inner ring is about 
19". There is also a nuclear ring with a semimajor axis ra- 
dius of about 4". The four bright spiral arms extend between 
20" and 40". The arms start from the inner ring, two from 
its major axis and the other two from the minor axis. Both 
the inner ring and the four spiral arms can be seen also in 
the NIR ff-band image. There is also a single outer spiral 
arm or a partial outer ring, that is not connected to the 
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Table 1. The main parameters for the three model families D, M 
and H. The H-band mass-luminosity ratios of the bulge (A//L)b 
and the disk (M/L)rj are given in solar units, assuming M h q = 
3.37 and the distance being 45 Mpc. Also listed are the parameters 
determining the halo rotation curve: velocity in infinity Voo in 
kms - 1 and core radius r c in arcseconds. 



Model family 


(A//L) B 


(M/L) u 


Voo 


r c 


D 


1.3 


1.3 






M 


1.2 


1.2 


152 


40.0 


H 


1.2 


0.54 


180 


20.0 



four-armed spiral l|Buta et al]ll99£l) . Fig. □ shows B- and 
//-band images (as observed in the sky and deprojected) of 
ESO 566-24. 

The rotation curve inferred from the velocity field (see 
Fig. |2l shows signs of probable non-circular motions: in the 
inner parts it rises steeply, then it has a short flat part at 
about 120 km s . Outside this area, it rises steeply again to 
the peak value, after which it quickly drops to the constant 
level of about 190 km s _1 . It is probable that the peculiar 
inner rotation curve is at least partly caused by non-circular 
motions induced by the bar, since axisymmetric rotation 
curves constructed from any realistic combinations of bulge, 
disk and halo cannot produce such a structure. Also the 
zero velocity curve of the velocity field implies the presence 
of non-circular motion. However, the effect is not as strong 
as in IC 4214 because in ESO 566—24, the orientation of the 
bar is almost perpendicular to the line of nodes, and hence 
little of the radial compone nt of the vel ocity vector projects 
to the line of sight (see also lButal lT988). 



5 MODELLING 

We have constructed three major sets of models: series D 
without a halo component, series M with a moderate halo 
and series H with a dominating halo (see Table 0. In all of 
these model series, we have varied several parameters such 
as the bar pattern speed and the disk thickness. 

The left side frames of Fig. |2] compare the rota- 
tion curves of models of series D (no halo), M (moderate 
halo) and H (dominating halo) with a similar disk thick- 
ness (hz/h r — 1/6). Also shown are the circular velocity 
points determined fro m the observed veloc ity field using the 
method presented bv I Warner et al.l (Il973l) . If we a dopt the 
definit ion of a maximal disk (or disk+bulge) by ISacketd 
(1997): a disk is maximal if its rotation speed at radius 
r2.2 = 2.2 h r (usually close to the peak of the disk rotation 
curve) is 75-95% of the total circular speed at this radius, 
then our moderate halo series M resides near the lower limit 
of maximal disks. Perhaps a better name would be a mini- 
mum halo model: only such an amount of halo was included 
that is needed to make the outer curve as flat as the observed 
rotation curve. 

The right side frames of Fig. show the correspond- 
ing "frequency diagrams", i.e. curves showing Q, Q ± k/2 
and ± ft/4, where £1 is the circular frequency and k is 
the epicyclic frequency. The inner parts of the models have 
very similar rotation curves, and thus with the same pat- 
tern speed, the inner resonance distances (ILRs and inner 



4/1) would be about the same. On the other hand, the outer 
resonance radii are different. For example, the distance be- 
tween the inner and outer 4/1-resonances is larger in models 
with a halo than in models without it. Throughout this ar- 
ticle we use km s _1 arcsec _1 as a unit of pattern speed. We 
denote models of each mass model series with the pattern 
speed: for example model D6.7 is a series D model with 
pattern speed Qt = 6.7 km s _1 arcsec _1 . At the adopted 
distance of 45 Mpc, 1 km s _1 arcsec _1 corresponds to about 
4.6 km s _1 kpc _1 and the length of a simulation is about 2.7 
Gyr (about 14 bar rotation periods in our best-fitting case) . 

5.1 The effect of model parameters 

5.1.1 The effect of disk thickness 

One of the main model parameters is the disk thickness. 
We have made simulations with a range of the disk scale 
height from 1/12 to 1/2 of the disk scale length, which should 
cover the plausible range of values for the disk thickness. 
In Fig. [3]we show the Fourier components of the tangential 
force relative to the axisymmetric radial force (total of bulge, 
disk and halo) as a function of radius for different mass 
models. 

As can be seen in Fig. [3] the m — 2 component has one 
maximum, at about 12", well inside the bar. The m — 4 
component has two clear maxima, one coinciding with the 
m = 2 maximum, and the other broader one in about the 
middle of the four-armed spiral structure. The m = 6 and 
m — 8 components have also two maxima, the inner coin- 
ciding with the m = 2 maximum, and the outer ones being 
inside the outer m = 4 maximum. 

As a measur e of bar strength, we use the relative tan- 
gential force Qt JCombes fc Sandersl[l98lh : 

Q T (r)=Fr x (r)/<F R (r)>, (6) 

where /""(r) is the maximum of tangential force at a given 
radius r, and < Fr(v) > is the azimuthally averaged radial 
force at the same radius. The maximum of Qr(r) over ra- 
di us, or Q;,, has been us ed as a single measure of bar st rength 
bv lButa fc Block! (1200 ill and lLaurikainen et alJ J2002I) . From 
Fig.|H]we can see that in our basic mass models, the extreme 
values of the Qt-parameter vary from 0.1 to 0.47 when going 
from the thickest disk models of series H to the thinnest disk 
models of series D. In each mass model series, the models 
with the thinnest disks have about twice as strong bars as 
models with thickest disks if measured by Qb. One can also 
see from Fig. |3] that the effect of disk thickness on Qt is 
much weaker in the region of the spiral arms. 

In Fig.|U we show the effect of disk thickness on the sim- 
ulated gas morphology. When comparing with the observed 
four-armed spiral of ESO 566-24, we see that in the D-series 
the fit is best when the disk is very thick (h z /h r — 1/2). 
In series M, the best fit is obtained with a thinner disk 
(both h z /h r — 1/12 and h z /h r — 1/6 give a fairly good 
fit). In series H, the spiral response to the potential is al- 
ways too weak. This can be explained by the weakness of 
the non-axisymmetric perturbation at the spiral region: it is 
less than half of that in series' M and D. We can also see 
that the size of the nuclear ring depends on the disk thick- 
ness: a thicker disk leads to a larger ring due to the reduced 
non-axisymmetric perturbation at its location, shifting the 
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Figure 2. Comparison between three models with different mass distributions (see TablelTl . The frames on the left side show the rotation 
curves, total rotation curve indicated as a continuous line, bulge as a dotted line, disk as a dashed line and the halo as a dash-dotted 
line. The asterisk symbols show circular velocities determined from the observed velocity field (see text). The frames on the right side 
show Q (thick continuous line), Q±k/2 (thin continuous lines) and Q±ft/4 (dotted lines) The dashed horizontal line shows the pattern 
speed at which the semimajor axis of the inner ring (19") coincides with the inner 4/1-resonance. The M/L of the bulge was chosen in 
each case so that the rotation amplitude at 20" corresponds with the observed one. 



ring closer to linear ILR. The axial ratio of the inner ring 
is not very sensitive to Qj, which is not surprising, because 
the location where the maximum of Qt is reached is well 
inside the ring. At the ring region, Qt varies less than its 
maximum value Qt when compared in models with different 
disk scale heights. 



5.1.2 The effect of pattern speed 

For each basic mass model, we ran a series of simulations 
with five different pattern speeds. The middle pattern speed 
of each series was chosen so that the semimajor axis radius of 
the inner ring coincides with the inner 4/1-resonance, and 
the highest pattern speed so that its corotation radius is 
1.1 times the bar radius. The other values of the pattern 



speed were chosen around the middle one with a step size 
1.2 km s _1 arcsec _1 . 

Fig-Elshows the effect of the pattern speed on gas mor- 
phology in series D. When the pattern speed rises, the nu- 
clear ring becomes smaller. With the highest pattern speed, 
Qb = 9.5, the inner ring disappears and the nuclear ring be- 
comes very elongated and parallel with respect to the bar. 
Four-armed spiral structure can be seen clearly in models 
D5.5 and D6.7. Both at higher and lower pattern speeds, 
the spiral tends to become two-armed. Several models in this 
and other mass model series have weak overdensities of gas 
particles in the bar area resembling leading offset dust lanes. 
These features are stronger in the early phases of simulations 
and more or less disappear during the simulations. These 
features move closer to the bar when the pattern speed is 
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Figure 3. The relative tangential force in models without a halo (top row), with a moderate halo (middle row) and with a strong halo 
(bottom row). The thick continuous line shows the maximum, whereas the m=2 component is indicated with a thin continuous line, 
m=4 component with a dotted line, m=6 with a dashed line and m=8 with a dot-dashed line. Note that owing to phase differences, the 
maximum is less than the sum of the amplitudes of the different components. 



increased, but in general they are too vague to be reliably 
compared with the observed dust lanes. 

Fig. HJ shows the effect of the pattern speed on the gas 
morphology in series M, which has a moderate halo com- 
ponent. The basic trend is quite similar as in the models 
without a halo, but the four-armed spiral appears in a larger 
pattern speed range (models M5.5 - M7.9), and is also more 
pronounced, the arms being longer and sharper. 

Fig. |7| shows the effect of the pattern speed on gas 
morphology in series H, with a dominating halo. The non- 
axisymmetric force is too weak to produce sharp spiral arms, 
and thus only weak spiral features appear, exhibiting a four- 
armed structure in models H5.5 - H7.9. 

5.1.3 Miscellaneous parametric dependencies 

We have also checked the effects of various other model pa- 
rameters: the number of Fourier components that were used 



to calculate the gravitational potential, the initial gas distri- 
bution, the bar turn-on speed, the collision frequency, and 
the coefficient of restitution. All these have only a minor ef- 
fect when compared to disk thickness and the pattern speed 
of the bar. 

The inclusion of Fourier component m — 4 and the 
higher ones makes the four-armed spiral more pronounced. 
When only m — and m = 2 components are included, 
the response is not a clear m = 4 spiral, but resembles 
model M6.7 with a thick disk. The other mass models fol- 
low similar trends. We have also made simulations where we 
include the odd components of the Fourier decomposition. 
This causes some asymmetry to the models, especially in the 
region where spiral arms emerge from the inner ring, but the 
effect is small. 

We found that the bar turn-on speed can affect the 
innermost morphology: when the bar was turned on too 
abruptly (obtaining full strength in one or less bar rota- 
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Figure 4. The effect of the disk thickness in models from the disk-dominated scries D (top row), series M with a moderate halo (middle 
row) and the halo-dominated series H (bottom row). Compare to Fig. |2]showing the strength of the non-axisymmetric force field in each 
case. The width of the frames is 120". The models are shown one bar period after the bar has reached its full strength. 



tion periods), some of the nuclear rings collapsed toward 
the centre or became very elongated and parallel to the bar. 
This is probably due to the high amount of intersecting or- 
bits. When the bar reached its full strength more gently (in 
2-4 bar rotation periods), these gas particles were slowly 
"guided" to orbits perpendicular with respect to the bar. 
However, when the pattern speed was high enough, even 
the gradual bar turn-on time could not save the nuclear 
rings, and the short-lived rings were destroyed before the 
bar reached its full strength. The reason for such behaviour 



is that the number of intersecting orbits increases when fib 
is increased. This was confirmed by simulations using non- 
colliding test particles. 

We studied the effect of collision frequency by mak- 
ing selected simulations with different gas particle sizes. 
If extreme values are ignored, the collision frequency af- 
fects mainly the sharpness of features and their formation 
timescale. However, the main morphological features are not 
strongly dependent on the collision frequency. We also stud- 
ied the effect of the coefficient of restitution a. Changing a 




Figure 5. Pattern speed series for models without a halo compared with a B-band image of ESO 566-24. The circles indicate radii of 
15, 30 and 45". The width of the frames is 120". The models are shown one bar period after the bar has reached its full strength. 



between our standard value, 0.0, and 0.5 has only a minor 
effect on the morphology. The greatest difference is in nu- 
clear ring: adopting higher values of a produces a smaller 
and thicker ring. 

5.2 A search for the best fit 

We use here the results of the previously presented simu- 
lation series to deduce the probable range of the essential 
model parameters. 

5.2.1 Morphology 

The top row of Fig. |H| shows the time evolution of 
morphologically-selected best-fitting models from the D- 
and M-series (in series H the spiral response is always very 
weak). The models without a halo, including the one shown, 
tend to have a problem with the positions of the four spiral 
arms: the arms are not equally spaced in azimuth. Instead 
of a clear m = 4 spiral, the structure more resembles two 
pairs of arms, or doubled arms. Another fault is the lack 
of "morphological longevity" in the models of series D: the 
spiral structure disperses in just a few bar rotation periods. 
The models with a moderate halo have more equally spaced 



and longer lasting spiral structure. A common trend in their 
long-term evolution is that when the spiral arms weaken 
owing to gas flow, the minor axis arms become relatively 
stronger than the major axis arms. This causes some models 
to become essentially two-armed in their later phases. Both 
model series produce slightly too pointed and elongated in- 
ner rings. Because the overall fit of the spiral morphology is 
better and due to its longevity, we prefer model M6.7 over 
model D6.7. 

We have tried to improve the pattern speed estimate 
by making several simulations around model M6.7. Judg- 
ing by the spiral morphology and the size of the inner 
ring, we derived £1), = 6.9 ± 0.5 km s _1 arcsec _1 , or about 
33.7 ± 2.3 km s _1 kpc _1 . Fig. H compares this morphologi- 
cal best-fitting model M6.9 with the B-band morphology of 
ESO 566-24. The resonance radii corresponding the epicycle 
approximation are indicated as circles. The nuclear ring is 
between two inner Lindblad resonances (1.3" and 9.4") and 
the inner ring is close to the inner 4/1-resonance (17.9"). The 
four-armed spiral structure is between the inner and outer 
4/1-resonances (38.7"). The OLR radius (46.5") is outside 
the spiral structure. When compared with the deprojected 
image of ESO 566-24, we can see that the bar ends well 
before the corotation resonance radius (28.5"), near the in- 
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Figure 6. Pattern speed series for models with a moderate halo. The circles indicate radii of 15, 30 and 45". The width of the frames 
is 120". The models are shown one bar period after the bar has reached its full strength. 



ner 4/1-resonance. Using the bar radius 18" ± 2", this gives 
rcn/rbar = 1.58. Combining the uncertainties in the deter- 
mination of rbar and Qt, we finally get rca/vtar = 1.6±0.3. 

The morphology of two "stellar" populations consisting 
of non-colliding test particles is also shown in Fig. |3] Both 
the low and medium velocity dispersion components (neither 
of these components should be expected to represent the self- 
gravitating stellar populations but more as a tracer of dy- 
namics) have a strong contribution of orbits resembling the 
closed £2-orbits (aligned perpendic ular with respect to the 
bar) of analytical bar potentials flContopoulos fc Grosbol| 
1989). Here, we call these orbits :r2-like orbits instead of 
X2-orbits for two reasons. First, they are orbits of individ- 
ual particles in our simulations, not exactly closed orbits 
found by iteration. Second, our potential differs from an- 
alytical bar potentials, because it is calculated from NIR- 
photometry. Thus, our bar is not straight (this causes the 
orientation of the inner orbits to be tilted with respect to 
the outer parts of the bar) and the contribution of the spiral 
arms is also included. Of course, the situation is somewhat 
similar in real galaxies, which do not accommodate idealis- 
tic bars. In addition to £2-like orbits, the presence of zi-like 
orbits (aligned parallel with respect to the bar) can be seen 
in both populations. It is evident that there is no support 



for a strong bar structure outside the inner 4/1-resonance. 
In the medium dispersion case, there is a low density feature 
with a minor-to- major axis ratio of about 0.5, that is aligned 
with the bar and reaches almost to corotation. However, it 
is more like an oval or a lens than a bar. 

In Fig. 1101 we compare the nuclear morphology of low 
dispersion collisionless populations in two different bar pat- 
tern speeds with the corresponding gas morphology and we 
also show some selected orbits of individual stellar parti- 
cles. In model M6.9, there is a "stellar" nuclear ring, which 
almost reaches the outer ILR radius of the epicycle approx- 
imation. The gaseous nuclear ring is considerably smaller. 
A possible reason for this difference is the crossing orbits, 
where colliding particles cannot stay. The outer X2-like or- 
bits are crossing both a considerable population of a;i-like 
orbits and other a;2-like orbits. Furthermore, the outermost 
orbits do not form (almost) closing loops, although they re- 
tain a perpendicular orientation with respect to the bar. The 
region suitable for non-crossing ^2-like orbits corresponds 
quite well with the size of the gaseous nuclear ring. In model 
M9.8 with higher pattern speed, the a;2-like orbits do not 
reach as close to the linear outer ILR radius. Their domain 
is also much narrower, which could explain why the gaseous 
nuclear ring has the orientation of rri-like orbits. 
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Figure 7. Pattern speed series for models with a strong halo. The circles indicate radii of 15, 30 and 45". The width of the frames is 
120". The models are shown one bar period after the bar has reached its full strength. 



We followed particle orbits in the region of the four- 
armed spiral and found a large number of particles librating 
aroun d Lagrangian points L4 and L5 (Bi nnev fc Tremainel 
119871) . which are close to the bar minor axis. However, the 
spiral structure cannot be explained by any single orbit type. 



Fig. 1111 shows the azimuthal variations of relative sur- 
face density at different radii both in observations and in 
our best-fitting model. We can see that although the m = 2 
component is dominating near the bar end, there are soon 
signs of an m = 4 pattern. In the midst of the spiral struc- 
ture, at about 35", all four arms are of about equal strength 
in the ff-band image, and the arm-interarm density contrast 
can be locally almost as high as 2.0 (about 0.7 magnitudes). 
Such high values in NIR-i mages have been fo und in only a 
few galaxies, such as M 51 (iRix fc Riekelll993D . Comparison 
between the gas component of model M6.9 and the B-band 
morphology of ESO 566-24 shows that the locations of in- 
tensity maxima are quite well fitted in the spiral arms. In 
the stellar component, the match is worse, which is not sur- 
prising since the model is not self-consistent. 



5.2.2 Kinematics 

In our previous modelling object, IC 4214, the fine details of 
the zero velocity contour of the velocity field were valuable 
in determining the bar strength and other model param- 
eters. Even though the non-axisymmetric perturbation in 
ESO 566-24 is considerably stronger than in IC 4214, the 
zero velocity contour is rather featureless. This is due to the 
unfavourable orientation of the bar with respect to the line 
of nodes. 

On the other hand, the "r otation curves" cons tructed 
using the method introduced bv I Warner et alJ lll973T) shows 
the effect of non-circular velocities. The clearest features 
affected by the non-axisymmetric perturbation are located 
inside 10", just inside the peak of Qt, the relative non- 
axisymmetric force. Furthermore, the flatness of the outer 
rotation curve clearly indicates that at least some amount 
of halo is needed (supposing constant M/L-ratio). 

Fig. 1121 shows comparison of selected rotation curves 
constructed from the modelled velocity field (data from 10 
time steps are summed) with the observed one. The observed 
velocity field does not have velocity information from the 
whole disk area, for example, the four-armed spiral struc- 
ture is incompletely covered. On the other hand, many of 
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Figure 8. The time evolution of two models from series D and M. The width of the frames is 120" and time is shown in simulation Gyr. 
The timesteps of the models are the same, but the times differ due to different mass scaling. 



our models include gas particles practically in the whole disk 
area. This is potentially dangerous when it comes to com- 
paring the observed and modelled kinematics. To diminish 
this problem we adopted a density threshold to the modelled 
velocity field. In some cases, this can have a notable effect 
on the inner part of the modelled rotation curve. 

Our morphologically fitted "standard" model, M6.9 
(frame a), gives a rather good fit with the position of the 
velocity jump near 10". The inner parts show also a simi- 
lar trend as the observed kinematics. However, the rotation 
amplitudes do not fit as well. Whereas the outer rotation 
amplitude is well fitted, the modelled velocity is too high 
inside 25". A model with a thicker disk (frame b) gives a 
better fit, but the rotation amplitude inside 6" is still too 
high, by about 40 km s -1 . This problem cannot be solved 
by changing the disk thickness or the pattern speed. A con- 
siderable increase of pattern speed destroys the two-stage 
structure of the rising part of the rotation curve, as is high- 
lighted by the kinematics of the model M9.8 (frame c). 

To get better kinematical fit, we tried changing the 
relative masses of bulge, disk and halo (the radial pro- 
files were not modified). In model B6.3 (frame d), where 
(M/L) B = 0.86, (M/L)t> = 1-02 and = 219 km s" 1 we 
could get the difference between the inner velocities as low as 
10 — 20 km s _1 . With the pattern speed 6.3 km s~ arcsec - , 
which sets the inner 4/1-resonance to approximately same 
distance as in model M6.9, the inner and nuclear rings have 
about the same size as in model M6.9. However, the four- 
armed spiral is not well fitted: instead of long spiral arms, we 
get several threaded arms. This is due to the increased halo 
contribution. If we further decrease the relative contribution 
of the bulge in order to reach the observed rotation ampli- 
tude, the nuclear ring changes its shape and orientation, 
by becoming dominated by ii-like orbits. Its morphology is 



then rather similar to the nuclear ring in model M9.8, and 
the central kinematics shows a peak instead of a drop. Thus, 
combining the morphological and kinematical comparison, 
model M6.9 is our best-fitting solution for ESO 566-24. 

The behaviour of the determined rotation curve inside 
10" is related to two factors: the orientation of the bar, i.e. 
its major axis is close to the kinematical minor axis, and to 
the orientation of x\- and a^Tike orbits with respect to the 
bar. In the case of pure elliptical orbits with major or minor 
axis parallel with the kinematical major axis, the resulting 
rotation curve is either below (orbit's major axis parallel 
with the kinematical major axis) or above (orbit's minor 
axis parallel w ith the kinem atical major axis) the local cir- 
cular velocity dTeubenlE )02'l. In barred galaxies, the shapes 
of orbits deviate from pure ellipses, but this approximation 
is still helpful to understand the general trends in the kine- 
matics of ESO 566-24. Now, the low velocities in the region 
inside 6" are caused by a^Tike orbits, whose major axis is 
perpendicular with respect to the bar major axis. The veloc- 
ity of these orbits near the bar minor axis (the major axis of 
the orbits), is much lower than the local circular velocity of 
the azimuthally averaged rotation curve. This is also close to 
the kinematical major axis, and thus the determined rota- 
tion curve in this region is below the axisymmetric rotation 
curve calculated from the mass distribution. 

Similarly, the velocity peak near 10" is caused by x\- 
like orbits, which are parallel to the bar. In this distance, 
they are near their minor axis, and the velocity is higher 
than the local circular velocity. The location of this feature 
depends on the extent of X2-like gas particle orbits. This is 
in part dependent on the pattern speed and the amplitude of 
the nonaxisymmetric perturbation. Also evolutionary effects 
are present: the jump can move inwards when the size of 
the nuclear ring decreases. In some simulations the jump 
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Figure 9. The resonance radii vs. morphology. The top left frame shows the resonance radii of model M6.9 plotted over the deprojected 
B-band image of ESO 566-24, in the top right frame they are plotted over the gas morphology of the model. The bottom frames show 
the morphologies of two populations of non-colliding test particles. In the low dispersion case, the velocity dispersion is about 1% of the 
local circular velocity, while in the medium dispersion case it is about 30% of the circular velocity. The circles drawn with a continuous 
line show the inner Lindblad resonance, the corotation resonance, and the outer Lindblad resonance, whereas the circles drawn with a 
dashed line shows the inner and outer 4/1-resonance radii. To increase resolution, the simulation images are made by summing particle 
positions of 10 different timesteps after the bar has reached its full strength. The width of the frames is 120". 



disappears when this region is depopulated of gas particles. 
The importance of £2-like orbits is further demonstrated by 
the kinematics of models which lack them: the plateau is 
replaced by a strong peak (M9.8 for example). 



6 DISCUSSION 

6.1 Spiral morphology in barred galaxies 

The origina l form of the Hubble classification of galax- 
ies (iHubbld Il926h divided spiral galaxies into subclasses 
a, b and c, depending not only the relative size of the 
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Figure 10. Stellar and gaseous nuclear rings vs. a;2-like orbits. The top row shows model M6.9, the left frame shows the low velocity 
stellar particles, the middle frame the gaseous nuclear ring, and the right frame selected stellar particle orbits. The bottom row shows 
the same for model M9.8. Outer and inner ILRs are indicated by circles drawn with a dashed line. The width of the frames is 20". 



bulge, but also on the pitch angle of the spiral arms and 
the degree of resoluti on in the arms. Sub s equent modifi- 
cations or extensions ijde Vaucou leurslll959t ISandagelll96lt 
ISandaee fc Bedkdll994l ; lButalll995l) . have partially removed 
the limitations of the original classification, e.g. by pro- 
viding description of ring and pseudoring morphologies. A 
dif ferent kind of spiral class i fication sch eme was developed 
bv lElmegreen fc Elmegreenl (|l982l ll987Fl . Simplifying their 
classification, most spiral galaxies can be divided into grand 
design, multiple-armed, and flocculent spiral galaxies. Grand 
design spirals have a global, usually two-armed spiral struc- 
ture. In multiple-armed galaxies, the spiral structure starts 
as two-armed, but the arms bifurcate later to several outer 
arms. Flocculent spirals have fragmented or chaotic spiral 
structure. 

A classic case of a barred spiral galaxy is NGC 1300, 
where a two-armed spiral structure starts from the ends of 
the bar. Quite often, the spiral arms emerge from an inner 
ring that surrounds the bar. There are also cases where the 
spiral arms start offset with respect to the ends of the bar. 
This can happen both in the leading and trailing quadrants 
of the bar (supposing that the spiral structure is trailing). 
Two-armed grand design structure is not the only possibility 
in barred galaxies. Different kinds of arm multiplicities also 
exist: arm bifurcation at large radius, doubled arms (one or 



both), and even genuine multiple-ar m structure . Flocculent 
ba rred galaxies are rare but known l)But,alll99fJ) . According 
to lElmegreen fc Elmegreenl 1)19951) , the length of the two- 
armed phase correlates with bar length: it reaches about 
twice the bar radius. 

Most work on galaxy morphology has been done using 
B-band plates. A major weakness of this is that the blue 
wavelength region is affected by dust extinction and the 
young stellar population, which usually constitutes only a 
minor part of the mass distribution. On the other hand, 
in the NIR most of the light comes from the old stellar 
population. It has been shown that bar s are more often de- 
tected in the NIR feskridge et al.l200Cil . and there are cases 
where the multiple-arm or flocculent spiral morphology of vi- 
sual images disappea rs and a two-armed grand design spiral 
emerges in the NIR llBlock fc Wainscoadll99ll: iBlock et alJ 
Il996t TT9991) . It has sometimes been claimed that this kin d 
of dualistic morpho logy is univer sal llBlock fc Puerarll9 99T 
However, lEskridge et alJ i2002f) . who used a large galaxy 
sample, demonstrated that dramatical differences between 
visual and NIR morphology are exceptions, and that in gen- 
eral galaxy morphology is quite similar in these two wave- 
length regions. 

The study bv lEskridge et alJ (l2002ft . based on the Ohio 
State University Bright Galaxy Survey (hereafter OSU- 
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Figure 11. The azimuthal variations of surface density at different radii (normalized to mean density at this distance) constructed from 
the Fourier decomposition of images. Each curve is calculated for the radius indicated by a line below it. The phase angle of the bar is 
zero. The two frames on the left side compare non-colliding test particles (low dispersion case) with the //-band image, and the frames 
on the right side compare gas particles with the B-band image. 



survey) , includes several barred galaxies where the multiple- 
arm structure is seen also in NIR images. The examples 
include all variations of arm multiplicities from arm bifur- 
cations to arm doublings and genuine multiple-armed spi- 
rals. For example, there are three-armed spiral galaxies, 
e.g. NGC 5054 (OSU-su rvey, another example is NGC 7137 
iGrosbol fc PatsisUl999tl . Although there are several four- 
armed galaxies in the OSU sample, none of them is similar 
to ESO 566-24. Most cases are arm bifurcations at large ra- 
dius or doubled arms, which are close to each o her. Even in 
cases where the extra arms do not form as bifurcations or 
doublings of an m — 2 spiral, they are usually weaker than 
the dominating two arms. 

ESO 566-24 with its very regular four-armed spiral 
structure is really exceptional. The four arms have roughly 
the same length, so they do not appear to have formed as 
a bifurcation of a usual two-armed structure. Furthermore, 
they are so far from each other (at 35", the phase differ- 



ence between the arms is about 90°), that arm doubling is 
clearly not the case. To crown the status of ESO 566-24 as 
a genuine m = 4 spiral, the arms have comparable strength 
also in the NIR. Thus, one cannot choose any pair of arms, 
say those starting near the bar ends, as two major arms. In 
the OSU-survey, the galaxy with the closest resemblance to 
ESO 566-24 is NGC 613. It has altogether five spiral arms, 
but the "extra" arms are not as strong and regular as in 
ESO 566-24. 

Some barred galaxies have features called plumes, i.e. 
short arcs besides the main spiral arms in the leading side 
with respec t to the bar ||Buta|l984j) . The best-known case is 
NGC 1433 feutall986l:lRvder et alll996ft . It is possible that 
plumes could be related to four-armed spiral structure, for 
example, they could be remnants of additional spiral arms. 
However, such features can form in gas dynamical simula - 
tions as result of doubled spiral arms fevrd et al]ll994D . 
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Figure 12. Comparison between the observed an d modelled kinemat ics. The continuous line shows the observed "rotation curve", 
derived from the velocity field using the method by War ner et alj Il973l) ■ The dotted line shows the axisymmetric rotation curve based 
on the mass model and asterisk symbols show the simulated rotation curve constructed in the same way as the observed one. 



It would be interesting to thoroughly compare ESO 566-24 
with NGC 1433. 

Although our best-fitting models can reproduce the 
four-armed spiral rather well, and give constraints to essen- 
tial model parameters, they do not directly unveil the origin 
of such structure. Instead, even in the best case, our models 
are snapshots of the current stage of the evolution of ESO 
566-24. It is possible that the location of resonances can give 
a clue of the origin of the four-armed spiral. Based on our 
models, it seems that four-armed structure is approximately 
confined between inner and outer 4/1-resonances. Thus the 
situation could be analogous to two-armed spirals, which are 
some times suggested to end near OLR fe.g. lElmegreen et alJ 
1989). Note however, that at least for non- barred galaxies, 
other solutions, which put the end of the spiral structure 
inside corotatio n or even the inner 4/1-r esonance, have also 
been suggested JPatsis et aUll994ll997h . 

6.2 Rings 

A connection bet ween rings and resonances w as perhaps 
first suggested bv lSchommer fc Sullivan! lll976t) . who iden- 
tified two rings in NGC 4736 with inner and outer Lindblad 
resonances. L ater, sticky particle simulations by ISchwarJ 
il98li Il984hl3) showed that gas can accumulate near the 



resonance radii. These and other studie s llCom bes fc Gerinl 
ll985tlBvrd et alJll994t iRautiainen fc SaloiboOOl) have iden- 
tified nuclear rings with orbits near the ILR, inner rings near 
the inner 4/1-resonance and outer rings near the OLR. Two 
differ ent observe d populations of outer rings and pseudor- 
ings |f3uta 1995) are nicely explained by orbits inside and 
outside OLR: inside OLR the major axes of the orbits be- 
longing to the dominating orbit family are oriented perpen- 
dicular to the bar, whereas outside OLR they are oriented 
parallel to the bar. 

In many studies, the nuclear rings in single bar poten- 
tials have been identified with a;2-orbits, which are perpen- 
dicular to the bar. When the non-axisymmetric perturba- 
tion is weak, and the deviations from circular motion are 
small, the outermost stable a;2-orbits reach the ILR (or the 
outer ILR when there are tw o inner Lindblad resonances) 
JContopoulos fc Grosbe)]ll989t) . However, when the strength 
of the perturbation is increased, the linear approximation is 
not valid, and the outermost extent of £2-orbits moves in- 
wards. In fact, some authors have redefined the ILR-radius 
so that it means the radius w i th th e outermost stable xi- 
orbits (Ivan Albada fc Sanders! Il983 lAthanassouIal Il992al : 
iMiwa fc Noguchilll998l) ! 

Recently. iRegan fc Teubenl (12003!) have studied the for- 
mation and evolution of nuclear rings using a hydrodynam- 
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ics code with an analytic bar potential. Their conclusion is 
that the nuclear rings are not related to inner Lindblad reso- 
nances, but instead they form as an interaction between gas 
streaming along 22- and ii-like streamlines. Furthermore, 
the nuclear rings in their simulations did not have a steady 
size but shrank throughout the simulation. Such behaviour 
can be seen also in our simulations, but the change in the size 
of the nuclear ring becomes very slow in the later phases of 
the simulations. This can be due to lack of interaction with 
particles on a;i-like orbits, which have become depopulated. 

We have tested the hypothesis by Regan and Teuben 
that the nuclear rings form as interaction of two different 
gas streamlines by making simulations where all the parti- 
cles initially resided in the domain of a^Tike orbits. A clear 
nuclear ring formed without interaction of two streamlines. 
In fact, a nuclear ring may form even without any gas: the re- 
distribution of non-colliding test particles by a rotating bar 
also forms a "concentration" of particles near the outermost 
T2-like orbits. Rings can form als o in self-gravitating stel- 
lar particles of N-body simulatio ns dRautiainen fc Salol200ol ; 
lAthanassoula fc Misiriotlsll2002l) . In ESO 566-24, the inner 
ring is clearly present in the H-band, and co rresponding 
examples exist also for outer and nuclear rings. iBvrd et alJ 
(1994) speculated that stellar rings could be remnants of 
star formation in gaseous rings, but it is also possible that 
a ring forms as a response of the old stellar population to a 
non-axisymmetric potential. 

If the rings are connected to orbits in the presence of 
a non-axisymmetric bar potential, then the s hapes of the 
rings should correlate with the bar strength. iButal (2002) 
compared the shapes of the inner rings in a few galaxies 
with the Q(,-parameter, and did not find a clear correlation. 
Galaxies with roughly similar Qb values can have inner rings 
of different shapes, ranging from virtually circular to highly 
elongated. This can be explained in several ways. First, Qb 
measures the peak strength of the bar, which often takes 
place well inside the inner ring, i.e. Qb may not be represen- 
tative of tangential force in the ring region. This effect can 
be seen in action in Figs. 5 - 7: in the higher pattern speed 
models, the inner ring is located closer to centre, in a region 
with higher Qt, which explains the more elongated shape 
(Qb is constant in each pattern speed series). Second, the 
Qb values were calculated from NIR photometry, and thus 
the possible contribution of dark matter was omitted. 



Tremaine- Weinberg method tend to be close to 1, although 
lAguerri et alJ (2003) found considerably higher values for 
galaxies NGC 1440 (1.6to.f) and NGC 3412 (1.5±g;|). Fur- 
thermore, the Tremaine- Weinberg method has been used al- 
most exclusively with SBO-galaxies. ESO 566-24 is of Hubble 
type SB(r)b, later than most galaxies whose pattern speed 
has been determined with the Tremaine- Weinberg method. 
lElmegreen fc ElmeereerJ dl985l) have suggested that early- 
type galaxies have longer bars than late-type galaxies in 
both the physical size and with respect to resonance loca- 
tions. Thus, it seems possible that the high value found in 
our study could somehow be related to Hubble type. 

A point in favour of higher pattern spee ds is the exis- 
tence o f dust lanes inside the bar. According to lAthanassoulal 
(1992b) the offset dust lanes form when the galaxy has an 
ILR and their shape corresponds to the observations when 
the corotation radius is 1.2 ±0.2 times the bar radius. How- 
ever, the dust lane argument for fast b a rs ma y not be as 
strong as it seems: iLindblad fc Krister] dl996f) found two 
possible pattern speeds for NGC 1300, giving rcR/rbar = 1.3 
and rcn/rbar = 2.4, both producing good dust lane mor- 
phology. 

N-body studies have shown that pattern speeds of 
the bar and spiral arms can be different, structures 
being either separate or in a non-linear mode cou- 
pling lISelTwood fc Sparkd Il988t iMasset fc Tagged Il997l : 
iRautiainen fc S alo 1999). Ra utiainen fc Said II2000I) demon- 
strated that rings can form also in self-gravitating simula- 
tions with two or more modes with different pattern speeds. 
Sometimes, the presence of two strong modes caused cyclic 
evolution in the shape and orientation of the ring. At least 
in the case of nuclear rings, this ca n be attributed to orbit 
loops (iMacieiewski fc Sparkdl2000h in a pote ntial with two 
rotating components jRmiriain^n" et alj l2002). So, in princi- 
ple it is it possible that ESO 566-24 could be modelled with 
two modes: the inner region having a higher pattern speed 
and the four-armed spiral a pattern speed in the range of our 
single pattern models. However, in this case the low value 
for the pattern speed is based not only on the modelling 
of the four-armed spiral structure, but also the central mor- 
phology is better modelled with low pattern speed: the inner 
and nuclear rings have approximately the right size. When 
the pattern speed is higher, both become too elongated or 
disappear altogether. 



6.3 Pattern speed dilemma 

Our best-fitting models have pattern speeds which place 
the corotation resonance at about 1.6 times the bar ra- 
dius. This is considerably larger than the often cited value 
rcR/r bar = 1.2± 0.2. For IC 4214 we found a value of 1.4±0.2 
dSalo et al Il999t) . thus just barely fitting within the fast bar 
regime although it has an early-type morphology. This raises 
a question if our modelling method is biased so that it gives 
too high values. Our preliminary results with NGC 4314 
suggest otherwise; for this galaxy we get rcn/rbar close to 

one. 

We extended a previous study bv lElmegreenl dl996T) by 
a brief search of the literature. We found pattern speed 
determinations of individual barred galaxies giving values 
ranging from below 1 to over 2, the average value being 
about 1.35 (median value was 1.3). The values based on the 



6.4 Halo contribution 

Early N-body simulations of barred galaxies dSerrwoodll98ll) 
exhibited deceleration of the bar rotation due to interaction 
with the outer disk. This made the corotation radius to move 
outwards, but it was compensated by the growth of the bar, 
i.e. the bar always ended close to corotation. However, in- 
teraction between bar and dark halo can introduc e bar slow- 
down that cannot be compensated by bar growth. IWemberg| 
( 1985) studied the dynamical friction between a bar and halo 
analytically and also by a "semi-resticted" N-body simula- 
tion (particles are not self-gravitating), and found that it 
causes the bar to lose most of its angula r mom entum in just 
a few bar rotations. iLittle fc Carlbergl lll99 J) studied bar 
slowdown by self-gravitating N-body models, where both the 
disk and halo were two-dimensional. The bar deceleration in 
their model was considerably lower, only by a factor of two 
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during 10 Gyr, half of which was due to interaction with 
the ou ter disk. On the other hand, iHernauist fc Weinberg! 
(1992) who used an analytical bar with a self-gravitating 
halo, found a slowdown rate corresponding to that found by 
Weinberg. 

At least seemingly confusing results have been obtained 
also with fully self-consistent N-body models o f bar-halo in- 
teraction. IDebattista fc Sellwoodl (ll99Sl l200Ch find that in 
models with centrally concentrated haloes, the bar slows 
down dramatically (rcR/rbar > 1-4), unless the halo has un- 
realistically high angular momentum in the same direction 
as the disk. IValenzuela fc Klvpinl l|2f)02h made simulations 
using adaptive grid refinement (the potential grid was sub- 
divided in high density regions) , which increased the gravity 
resolution. They found that the bar pattern speed was al- 
most constant and that the stellar disk lost only 5 - 10% of 
its angular momentum to the halo. Furthermore, they could 
reproduce the dramatic bar slo wdown when the gri d reso- 
lution was decreased. Recently. lAthanassoulal i2003T) found 
that when the halo velocity dispersion is increased, the bar 
slowdown rate also decreases. 

The preceding discussion indicates that it is question- 
able whether bar pattern speed can be used as an argu- 
ment for high or low halo-to-disk mass ratios. If we com- 
pare our best-fitting mass m odel with the simulations of 
IDebattista fc Sellwoodl 1^998) with a non-rotating halo, and 
adopt their parameter r\ = (vdisk/vhaio) 2 at the maximum 
of disk rotation curve, we find rj — 3.65. In their Fig. 2, 
this would correspond to rca/rbar = 1.3 — 1.4 at the fi- 
nal equilibrium state. On the other hand, if we include 
bulg e to Vhaio, we get y = 1 . 39, co r responding r C R/r bar > 
2 in IDebattista fc Sellwoodl dl99gl) . IDebattista fc Sellwoodl 
(2000) found that a halo which rotates in the same sense 
as the bar can induce lower bar slowdown rate than non- 
rotating or retrograde halos. When taking into account that 
bulges rotate, our res ult real Thar ~ 1.6 is not in d i sagree - 
ment with results of IDebattista fc Sellwoodl il998L 12000). 
However, t he situation is the same if w e compare with sim- 
ulations bv lValenzuela fc Klvpinl ((2002) : in their models the 
corotation resonance radius is 1.2 to 1.7 times the bar radius, 
and our value is near the upper limit of their range. 

In our best fitting model (M6.9), the dark halo con- 
tributes less than 20% of the mass inside T2.2- The rota- 
tion curve is dominated by luminous matter even inside the 
whole disk region. Thus, our results with ESO 566-24 are 
in accordance with estim ates of disk mass co ntribution by 
dynamical modelling (e.g. IWeiner et al]l200ll) . but disagree 
with cold dark matter (CDM) cosmolo gical N-body simula - 
tions, which produce cuspy halos (e.g. iNavarro et alj|l996l) 



7 CONCLUSIONS 

We have constructed dynamical models for the four-armed 
barred spiral galaxy ESO 566-24. The mass distribution of 
the different components is based on near-IR photometry 
with the exception of a possible dark halo, whose amount 
was chosen to fit the rotation curve at large radii. We have 
been able to construct a model which reproduces the ob- 
served four-armed morphology. Also the main kinematic 
characteristics are seen in our models. The main conclusions 
are as follows. 



1. Two major factors affecting the simulated morphol- 
ogy and kinematics are the strength of the non-axisymmetric 
perturbation and the pattern speed of the bar. The former 
depends on the mass and the thickness of the disk. 

2. Both the kinematical and morphological fits are bet- 
ter when a moderate halo component is included. On the 
other hand, a dominating halo can be ruled out, because 
then the gas response to the disk gravitational potential is 
too weak. Thus, the contribution of the luminous matter 
dominates the rotation curve in the whole optical disk. 

3. The four-armed spiral could be produced between 
the inner and outer 4/1-resonances. The inner ring is close 
to the inner 4/1-resonance and the nuclear ring is between 
two inner Lindblad resonances. The size of the nuclear ring 
seems to correlate with the size of the region suitable for 
non-crossing 12-like orbits. 

4. If the bar has the same pattern speed as the four- 
armed spiral, then it rotates rather slowly; the corotation 
radius is about 1.6 ± 0.3 times the bar radius (including un- 
certainties both in the bar radius and pattern speed) , which 
is higher than the often cited value of 1.2 ± 0.2. The exis- 
tence of several other determined high values for this ratio 
suggests that there is more scatter in the rca/rbar ratio than 
is sometimes stated. It should also be noted that most bar 
pattern speed determinations have been made for galaxies 
of earlier Hubble type than ESO 566-24. 

5. The rather high value of rcn/rbar is difficult to be dis- 
puted by pattern speed multiplicity. In addition to the spiral 
structure, the nuclear and inner rings can be reproduced in 
the right scale with a single pattern speed. Adopting con- 
siderably higher pattern speeds, which would be consistent 
with the bar ending near corotation, produce significantly 
worse-fitting ring morphology. Also, the kinematics of the 
central parts deviate more from the observed velocities. 

6. The relatively slow bar rotation rate can be due to 
interaction between the bar and the spheroidal component 
(halo and bulge) . The exceptional m = 4 spiral morphology 
can be related to the low bar rotation rate and the impor- 
tance of the spheroidal component. 

7. The curious shape of the inner rotation curve can be 
explained by combined effects of xi- and a^-like orbits. 
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